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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

SHINY  CONTACT/DISTRIBUTED  FEEDBACK  SEMICONDUCTOR  LASERS 

By 

HORNG-JYE  LUO 
August  1993 

Chairman:  Peter  Zory 

Major  Department:  Electrical  Engineering 

Shiny  contact  diode  lasers  have  been  demonstrated  for 
the  first  time  in  the  1 . 3(lm  InGaAsP/InP  material  system. 
Power  levels  greater  than  one  watt  pulsed  have  been  achieved 
at  room  temperature.  Using  measured  parameters,  improvements 
in  continuous  wave  performance  levels  relative  to 
conventional  devices  are  predicted.  The  advantages  of 
combining  the  shiny  contact  with  the  distributed  feedback 
(DFB)  laser  structure  are  described.  Ray  optics  and  Floquet- 
Bloch  formalisms  are  developed  to  model  the  lasing  properties 
of  these  shiny  contact  DFB  lasers  as  well  as  conventional  DFB 
lasers.  It  is  shown  that  the  ray  optics  and  Floquet-Bloch 
formalisms  are  superior  to  the  standard  overlap  integral  and 
transfer  matrix  techniques  in  strong  coupling  configurations. 
The  ray  optics  formalism  is  extended  to  deep  groove  gratings, 
of  potential  use  in  fabricating  low  threshold  current, 
grating  surface  emitting  laser  arrays. 
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CHAPTER  I 
INTRODUCTION 


Semiconductor  diode  lasers,  due  to  their  small  size,  are 
extremely  "hot"  devices.  For  example,  the  intensity  at  the 
output  facet  of  a lmW  diode  laser  is  comparable  to  that  at 
the  surface  of  the  sun  [Zor92] . Therefore,  heat  dissipation 
is  a critical  issue  in  the  design  of  continuous  wave  (CW) 
diode  lasers  for  high  power  or  high  temperature  applications 
such  as  analog  communications,  optical  recording,  tunable 
spectroscopic  systems  and  pump  sources  for  solid  state 
lasers . 

The  maximum  CW  output  power  of  a diode  laser  is  limited 
by  two  factors:  catastrophic  optical  damage  [Hen79]  and 
thermal  roll-over  [Tod86] . The  performance  of  long 
wavelength  diode  lasers  such  as  InGaAsP/InP  is  usually  more 
sensitive  to  the  temperature  variation  than  that  of  short 
wavelength  lasers  such  as  AlGaAs / GaAs . When  these  long 
wavelength  lasers  are  operated  at  high  power  levels,  the 
internal  heating  raises  the  temperature  of  the  active  layer, 
causing  the  threshold  current  to  increase  and  efficiency  to 
decrease,  which  further  elevates  the  temperature.  Hence  the 
maximum  CW  output  power  is  usually  limited  by  thermal  roll- 
over. Since  the  temperature  sensitivity  for  double- 
heterostructure (DH)  lasers  is  basically  a material  parameter 
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and  cannot  be  changed  for  a given  material  system,  one 
possible  solution  to  increase  the  output  limit  of  thermal 
roll-over  is  to  lower  the  thermal  impedance  between  the  heat 
source  (active  layer)  and  the  heat  sink  (copper  block) . This 
can  be  done  using  the  shiny  contact  laser  structure  [Zor75] 
in  which  the  thermal  impedance  is  minimized  by  the  use  of  a 
thin  cladding  layer.  In  this  work,  the  first  operation  of 
such  devices  in  the  1 . 3|lm  InGaAsP/InP  material  system  is 
demonstrated  and  power  levels  greater  than  1W  pulsed  achieved 
at  room  temperature.  This  is  among  the  highest  power  levels 
ever  reported  for  single  element  InGaAsP/InP  DH  lasers  and  is 
comparable  to  those  produced  by  high  power  separate 
confinement  heterostructure  (SCH)  single  quantum  well  (SQW) 
InGaAsP  lasers  operating  at  1 . 3|lm  [Zor93]. 

By  combining  the  shiny  contact  with  the  distributed 
feedback  ( DFB ) laser  structure,  more  advantages  can  be 
expected  such  as  higher  device  reliability,  better  single 
frequency  yield  and  lower  fabrication  cost  compared  with 
conventional  type  DFB  lasers  [Luo90] . The  shiny  contact  DFB 
laser  is  different  from  conventional  DFB  lasers  in  that  the 
diffraction  grating  is  a corrugated  metal/semiconductor 
interface  as  opposed  to  a corrugated 
semiconductor/semiconductor  interface.  In  conventional  DFB 
lasers,  the  dielectric  grating  is  only  about  0 . l|lm  from  the 
active  layer,  so  the  defect  states  on  the  corrugated  surface 
may  propagate  into  the  active  layer  during  the  high 
temperature  epitaxial  regrowth  after  grating  fabrication.  In 
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the  shiny  contact  DFB  laser  configuration,  the  diffraction 
grating  is  placed  about  an  optical  wavelength  (0.3~0.5|lm) 
from  the  active  layer  and  the  ohmic  contact  is  made  directly 
on  top  of  the  grating.  Thus  the  fabrication  cost  is  reduced 
and  the  reliability  issue  is  avoided  since  no  epitaxial 
regrowth  is  required. 

In  modeling  laser  action  in  DFB  devices,  one  normally 
utilizes  coupled  mode  theory  [Kog72] . In  this  theory,  the 
key  laser  parameter  is  the  backward  coupling  coefficient  K 

since  its  magnitude  determines  such  parameters  as  frequency 
selectivity  and  threshold  mode  gain.  The  standard  technique 
to  determine  K in  a DFB  laser  is  to  evaluate  the  coupled  mode 

integral  which  describes  the  overlap  of  the  unperturbed  mode 
with  the  grating  perturbation  [Str75]  . It  has  been  shown 
that  this  technique  shows  an  inherent  problem  as  to  how  to 
choose  the  unperturbed  waveguide  geometry  when  the 
perturbation  is  strong  [Str77]  [Cor88].  For  example,  in 
semiconductor/metal  or  semiconductor/air  gratings  where  the 
refractive  index  difference  is  large,  the  standard  approach 
gives  results  strongly  dependent  on  the  choice  of  the 
unperturbed  waveguide  [Sha89].  Another  widely  used  method, 
the  effective  index  transfer  matrix  method  [Han92],  shows 
difficulties  in  considering  transverse  modes  and  ambiguities 
in  the  transformation  for  highly  asymmetrical  corrugated 
waveguides . 

As  a result,  ray  optics  and  Floquet-Bloch  formalisms  are 
developed  to  model  the  lasing  properties  of  the  shiny  contact 


4 


DFB  structure  where  standard  techniques  fail.  The  ray  optics 
formalism  is  also  used  to  derive  closed  form  expressions  for 
the  coupling  coefficient  in  conventional  type  DFB  structures 
with  dielectric  gratings  [Luo93].  In  such  weak  coupling  DFB 
configurations,  the  K values  obtained  from  the  ray  optics 

expression,  transfer  matrix  and  overlap  integral  approaches 
show  excellent  agreement.  The  ray  optics  formalism  is 
superior  to  standard  techniques  because  it  provides  device 
designers  with  closed-form  expressions  in  both  strong  and 
weak  coupling  DFB  configurations,  and  it  can  be  extended  to 
deep  groove  gratings. 

Grating  surface  emitting  (GSE)  lasers  can  be  integrated 
into  two  dimensional  arrays  and  are  widely  used  in  high  power 
applications  such  as  pump  sources  for  solid  state  lasers.  In 
the  conventional  design  with  shallow  groove  gratings, 
however,  the  threshold  current  is  high  due  to  the  presence  of 
the  radiation  loss.  Therefore  GSE  lasers  have  not  been 
utilized  in  low  power  applications  such  as  chip-to-chip 
communications  [Tak92].  Using  the  extended  ray  optics 
formalism,  the  first  low  threshold  design  of  GSE  lasers  is 
proposed  based  on  deep  groove  shiny  contact  gratings . These 
structures  could  be  competitive  with  vertical  cavity  surface 
emitting  lasers  (VCSEL)  [Jew91]  . 

This  dissertation  is  organized  as  follows.  In  the 
second  chapter,  the  fabrication  and  operation  of  1 . 3|lm 
InGaAsP/InP  shiny  contact  DH  lasers  are  discussed,  and  the 
experimental  data  successfully  explained  by  the  theory. 
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Using  the  parameter  values  extracted  from  the  measurements, 
the  improvements  in  high  power  CW  performance  levels  relative 
to  conventional  devices  are  predicted. 

Chapter  III  discusses  the  application  of  the  Floquet- 
Bloch  formalism  to  the  shiny  contact  DFB  laser.  The  general 
formalism  is  truncated  by  including  only  the  minimum  number 
of  terms  and  applied  to  several  special  configurations. 

In  Chapter  IV,  closed-form  expressions  for  the  DFB 
coupling  coefficient  are  derived  using  the  ray  optics 
formalism  in  shiny  contact  and  conventional  DFB  laser 
configurations.  Numerical  results  using  the  Floquet-Bloch, 
ray  optics  and  standard  techniques  are  compared  in  various 
DFB  structures . 

Then  the  deep  groove  effects  are  investigated  in  Chapter 
V.  The  ray  optics  formalism  is  extended  to  design  the  deep 
groove  shiny  contact  GSE  laser. 

Finally,  Chapter  VI  contains  conclusions  and 


recommendations  for  future  research. 


CHAPTER  II 

HIGH  POWER  InGaAsP/InP  SHINY  CONTACT  DIODE  LASERS 


It  is  well  known  that  the  experimentally  observed 
threshold  current  Ith  for  double  heterostructure  (DH)  diode 
lasers  has  the  following  temperature  dependence  [Cas78] : 

Ith  ~ exp (T/T0)  (2.1) 

where  T is  the  temperature  of  the  active  layer  of  the  laser 
and  To  the  characteristic  temperature.  Generally  speaking, 
long  wavelength  lasers  have  a low  value  of  Tq . For  example, 
To=170K  for  AlGaAs/GaAs  DH  lasers  and  Tq^SSK  for  InGaAsP/InP 
DH  lasers  around  room  temperature  [Agr86] . This  indicates 
that  the  performance  of  InGaAsP  lasers  is  much  more  sensitive 
to  temperature  than  that  of  AlGaAs  lasers.  To  explain  the 
stronger  temperature  dependence  for  InGaAsP  lasers,  various 
mechanisms  have  been  proposed  and  investigated,  such  as  Auger 
recombination  [Sug82],  leakage  current  [Yan81],  free  carrier 
absorption  [Asa81]  and  intervalence  band  absorption  [Asa83]. 
No  matter  which  mechanism  is  dominant  and  responsible  for  the 
lower  value  of  To,  this  phenomenon  causes  thermal  roll-over 
to  limit  the  highest  CW  power  achievable  for  InGaAsP/InP 
diode  lasers,  in  contrast  to  AlGaAs/GaAs  lasers  where 
catastrophic  optical  damage  is  usually  the  limiting  factor. 
This  is  due  to  the  elevated  temperature  in  the  active  layer 
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under  high  power  CW  operation  since  not  all  of  the  input 
electrical  power  is  emitted  as  light  but  a large  fraction  of 
it  is  converted  into  heat  . When  the  temperature  of  the 
active  layer  is  increased,  the  threshold  current  increases 
and  efficiency  decreases,  which  further  raises  the 
temperature.  Hence  the  thermal  roll-over  sets  a limit  for 
InGaAsP/InP  lasers  in  high  power  or  high  temperature 
applications . 

Since  the  characteristic  temperature  To  for  DH  lasers  is 
basically  a material  parameter  and  can  not  be  changed  for  a 
given  material  system,  one  possible  way  to  improve  the 
thermal  roll-over  limit  is  to  lower  the  thermal  impedance 
between  the  heat  source  (active  layer)  and  the  heat  sink.  In 
conventional  diode  lasers,  the  thickness  of  the  cladding 
layer  is  usually  1 . 5(lm  and  is  far  more  than  enough  [Beh90]  . 
In  this  chapter,  we  propose  and  demonstrate  a shiny  contact 
laser  structure  where  the  thermal  impedance  is  minimized  with 
the  use  of  a thin  cladding  layer.  Reduction  of  the  thickness 
of  the  top  clad  offers  several  advantages.  First  of  all,  the 
heat  generated  in  the  active  layer  is  more  efficiently 
removed  because  the  thermal  impedance  is  minimized. 
Consequently,  devices  of  this  type  can  produce  higher  power, 
operate  at  higher  temperature,  and  have  better  reliability 
and  durability.  Secondly,  a thin  cladding  layer  requires 
less  epitaxial  growth  time  and  source  material.  This  is 
especially  attractive  for  the  molecular  beam  epitaxy  (MBE ) . 
Other  advantages  become  obvious  when  the  shiny  contact  is 
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combined  with  the  distributed  feedback  (DFB)  laser  structure, 
and  are  presented  in  the  next  chapter. 

The  main  problem  of  such  laser  structures  is  the 
fabrication  of  the  shiny  contact,  which  must  have  both  low 
electrical  resistance  and  low  optical  loss.  We  show  that 
this  problem  can  be  solved  and  both  theoretical  and 
experimental  results  of  1 . 3(lm  InGaAsP/InP  shiny  contact  diode 
lasers  are  given.  Power  levels  greater  than  one  watt  pulsed 
have  been  achieved  at  room  temperature.  To  the  best  of  our 
knowledge,  this  is  the  first  demonstration  of  such  devices  in 
1 . 3|lm  InGaAsP/InP.  The  power  levels  we  have  achieved  are 
comparable  to  those  obtained  from  quantum  well  InGaAsP/InP 
lasers  [Zor93].  Section  2.1  discusses  the  theoretical  design 
of  the  shiny  contact  laser  structure,  followed  by  discussions 
of  the  fabrication  procedures  and  experimental  results. 
Finally  the  potential  impact  of  such  laser  structures  is 
predicted . 

2 . 1 Theoretical  Analysis  of  Shiny  Contact  Lasers 

Figure  2-1  shows  the  DH  type  InGaAsP/InP  laser  structure 
and  relevant  parameters  used  throughout  the  theoretical 
analysis . The  complex  refractive  index  of  gold  is  obtained 
from  [Pal85],  and  the  extinction  index  of  the  ternary  InGaAs 
is  derived  from  an  estimated  absorption  of  6700/cm  at  lasing 
wavelength  A.=  1.3[lm.  The  thickness  of  the  contact  layer  is 

the  minimum  value  required  for  a good  ohmic  contact  as 
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determined  by  the  depletion  width  of  the  metal-semiconductor 
junction.  Gold  is  the  natural  candidate  for  the  contact 
metal  due  to  its  superior  electrical  and  optical  properties. 
The  p-clad  thickness  t is  the  main  variable  whose  value  is  to 
be  optimized. 

To  investigate  the  effects  on  the  laser  performance  when 
t is  reduced,  the  two  key  factors  are  the  optical  mode 
confinement  factor  T and  the  optical  mode  loss  am  due  to  the 

metalizat ion  (introduced  by  the  contact  layer  and  the  contact 
metal)  . These  two  factors  directly  affect  the  threshold 
current  and  efficiency  of  a laser.  In  general  one  wants  to 
minimize  the  loss  but  maximize  the  confinement. 

The  optical  mode  loss  due  to  the  metalization  can  be 
calculated  by  solving  the  transcendental  eigenvalue  equation 
for  the  five  layer  waveguide  structure  [Tam90]: 

..  . . (P0P2+P12)  sinh (pih)  + pi (p2+po)  cosh (pih) 

exp(2p2t)  o 

(P0P2~Pi  ) sinh (pih)  + pi (p2-po)  cosh (pih) 

(P2P4-P32)  Sinh (p3d)  + P3 (P2~P4) cosh (p3d) 

(P2P4+P32)  sinh (p3d)  + P3 (P2+P4) cosh (p3d) 

where 

Pi  = VPo2  - ni2k2  i=0,  1, 2, 3,  4 

k = 27C/X,  X:  wavelength  in  vacuum 

Po  = neff  k nef f : effective  waveguide  index 

The  loss  due  to  the  free  carrier  absorption  in  the  metal  and 
the  band-to-band  absorption  in  the  contact  layer  are  included 
in  the  complex  refractive  indices,  and  therefore  the  optical 
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mode  loss  am  of  the  waveguide  is  directly  related  to  the 
imaginary  part  of  the  propagation  constant  pQ, 

(Xm  = 2lm(P0)  (2.3) 

In  Figure  2-2  the  calculated  results  for  am  are  plotted 

as  a function  of  the  p-clad  thickness  t.  It  is  observed  that 
(Xm  increases  monotonically  as  t decreases.  However,  since 

the  internal  optical  loss  for  1 . 3|lm  InGaAsP  DH  lasers  is 
typically  around  30cm_1,  the  extra  loss  introduced  by  the 
metalization  can  be  neglected  for  t>0.4|0.m. 

The  optical  mode  confinement  factor  T is  defined  as  the 
proportion  of  the  optical  mode  overlapping  with  the  active 
layer  where  optical  gain  is  present: 
d 

J | E ( x ) | 2 dx 

m 0 

r * — (2.4) 

Jl E (x)  1 2 dx 

— oo 

where  E (x)  is  the  electric  field  distribution  for  the 
transverse  mode  of  interest.  One  can  solve  E(x)  then 
numerically  integrate  eq.(2.4)  to  obtain  T.  Alternatively, 

one  can  put  a small  material  gain  Ag  in  the  active  layer  and 
solve  for  the  complex  propagation  constant  of  the  waveguide, 
then  r is  the  relative  increase  in  the  mode  gain  AG  per  unit 

Ag . The  latter  technique  is  utilized  to  obtain  the  results 
shown  in  Figure  2-3.  It  is  seen  that  as  t is  reduced,  T 

increases  slowly  in  the  beginning  then  goes  down  suddenly. 
There  is  a maximum  value  for  T occurring  at  t«0.2p.m.  This 

phenomenon  can  be  attributed  to  the  competition  between  the 
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effects  of  the  shiny  contact  and  the  twin  guide  formed  by  the 
active  and  contact  layers.  When  t is  reduced,  the  metal 
tends  to  push  the  optical  mode  downward  while  the  formation 
of  the  twin  guide  tends  to  trap  the  optical  power.  The 
latter  effect  dominates  for  0 . 5(lm>t>0 . 2(lm . As  t is  smaller 
than  0.2|lm,  however,  the  effect  of  the  shiny  metal  becomes 
dominant  and  the  optical  mode  is  pushed  away  so  T drops 

sharply . 

The  threshold  mode  gain  Gth  and  the  threshold  material 
gain  gth  for  a laser  are 

1 1 

r gth  = Gth  = OCi  + Cim  + ^ In  (^) 


= dt  + l ln(|)  (2.5) 

where  (Xi  = 30cm-1  is  the  internal  optical  loss  due  to  free 
carrier  absorption,  at  the  total  optical  mode  loss,  L=250|Im 
the  cavity  length  and  R=0.32  the  facet  reflectivity.  The 
carrier  density  at  threshold  Nth  can  be  empirically 
approximated  as  the  following  function  of  gth  [Dut82]  : 

Nth  = 1.103X1018  + 6 . 445xl015-gth  - 3 . 856xl013-gth2 
+ 1 .VGSxlO^-gth3  - 4.15  6x1 0 8-gth4  + 3 . 993xl05-gth5  (2.6) 

where  Nth  is  in  1/cm3  and  gth  in  1/cm.  The  radiative  and 
nonradiative  current  density  can  be  written  as  [Dut82] 


Jr  ~ p ( 5th  + P jtr ) 

3 

Jnr  — e d C Nth 

where  P=0 . 0783cm-|lm/A  and  jtr=3000A/cm2-|lm  are 
parameters  and  C=3 . 0xl0“29cm6/s  the  Auger 


(2.7) 

(2.8) 

calculated  gain 
coefficient  for 
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InGaAsP  [Agr86].  The  threshold  current  density  Jth  and  the 
differential  quantum  efficiency  are  then 

Jth  = ( Jr  + Jnr)  = t|i  T|r  Jr  = TlSp  Jr  (2 . 9) 

l ln<£> 

Tld  = TU  ^ (2.10) 

where  TJi  is  the  injection  efficiency  defined  as  the  fraction 

of  carriers  injected  into  and  recombined  in  the  active  layer, 
TJr  the  radiative  efficiency  (or  internal  quantum  efficiency) 

defined  as  the  fraction  of  injected  carriers  recombined  as 
photons,  and  TlSp  the  spontaneous  efficiency  - the  product  of 
T)i  and  T|r.  These  parameters  can  be  extracted  from  the 
experimental  data  (see  Section  2.3)  and  we  assume  Tii  = l in 

this  section. 

The  calculated  threshold  current  density  and 

differential  quantum  efficiency  are  presented  in  Figures  2-4 
and  2-5  respectively.  It  is  interesting  to  note  from  Figure 
2-4  that  as  t decreases,  Jth  actually  goes  down  a little  bit 
before  it  rises  up  sharply.  This  is  due  to  the  fact  that  the 
increase  in  the  optical  mode  loss  is  more  than  offset  by  the 
improvement  in  the  confinement  factor  as  one  reduces  the  p- 
clad  thickness.  Since  this  effect  is  very  subtle,  one  can 
generally  assume  that  Jth  is  constant  at  t>0.4|lm.  The 

differential  quantum  efficiency  decreases  monotonically  as  t 
decreases  but  does  not  change  much  until  t<0.5[lm.  The 

reduction  of  Tjd  simply  comes  from  the  increase  in  at. 
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From  the  above  theoretical  analysis,  we  see  that  the  p- 
clad  thickness  in  conventional  diode  laser  is  much  larger 
than  needed.  In  principle  one  can  improve  the  high  power  CW 
performance  of  InGaAsP  DH  lasers  by  using  a thin  p-clad  of 
about  0 . 5flm  without  degrading  such  properties  as  threshold 
current  density  and  differential  quantum  efficiency.  In  the 
next  section  the  fabrication  procedures  of  such  shiny  contact 
lasers  in  InGaAsP/InP  are  discussed. 

2 . 2 Fabrication  of  Shiny  Contact  Lasers 

Wide-stripe  gain-guided  InGaAsP/InP  shiny  contact  lasers 
were  fabricated  from  Wafer  QA21125  grown  by  Lasertron  using 
liquid  phase  epitaxy  (LPE) . A schematic  diagram  for  the  DH 
laser  structure  is  shown  in  Figure  2-6  where  all  relevant 
parameter  values  are  given.  The  processing  techniques  are 
described  as  follows. 

The  SiC>2  coated  wafer  was  originally  about  400|lm  thick 
and  was  lapped  down  to  about  100|im  in  thickness  using  9|lm 
AI2O3  powder.  During  the  lapping  process,  the  wafer  were 
surrounded  by  3 GaAs  dummy  stacks,  each  of  which  contained  2 
GaAs  pieces  adhered  together  by  photoresist.  Since  InP  based 
materials  are  softer  than  GaAs,  a lower  lapping  speed  ~10rpm 
was  used. 

The  wafer  was  then  cleaned  in  heated  TCA,  acetone, 
methanol  and  IPA,  rinsed  in  DI  water  and  blown  dry  using  a 
nitrogen  gun.  After  dehydrating  at  110°C  for  30  minutes,  the 
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wafer  was  placed  on  a cover  glass  which  had  been  spun  coated 
with  photoresist  1350J  at  2000rpm  for  15  seconds.  The  glass- 
wafer  was  baked  at  90°C  for  10  minutes  so  that  the  wafer  would 

not  drop  during  processing.  Then  photoresist  1350J  was  spun 
coated  on  SiC>2  of  the  wafer  at  5000rpm  for  30  seconds,  and  a 
soft  bake  was  performed  at  90°C  for  20  minutes.  After  an 
exposure  of  10  seconds  with  the  mask  on  the  front  (SiC>2  side) 
and  1 minute  without  the  mask  on  the  back  (through  the  cover 
glass),  the  glass-wafer  was  developed  at  1:7.5  diluted 
developer  606  for  12  seconds  and  the  wafer  was  separated  from 
the  glass . The  wafer  was  rinsed,  blown  dry  and  hard  baked  at 
110°C  for  20  minutes  before  chemically  etched  in  BOE(6:l)  for 

2 minutes  and  45  seconds  to  open  a 100|lm  window  in  the  SiC>2 
layer . 

Immediately  following  another  rinsing  and  blowing,  the 
epitaxial  side  was  e-beam  evaporated  with  1500A  gold  at 
3A/sec  with  liquid  nitrogen  in  the  diffusion  pump.  A lift- 
off was  then  done  in  hot  acetone  so  that  only  the  gold  on  top 
of  the  window  (in  direct  contact  with  the  contact  layer)  was 
left.  After  a regular  cleaning  procedure,  the  whole  oxide 
side  was  coated  with  1000A  nickel  and  500A  gold.  Finally  the 
substrate  side  was  coated  with  2000A  gold. 

The  wafer  was  cleaved  into  bars  of  cavity  length  ranging 
from  125jlm  to  500|lm  using  our  standard  cleaving  techniques: 
scribe,  roll  and  break.  The  cleaved  bars  were  probe  tested 
before  they  were  diced  into  laser  chips.  Then  some  of  the 
chips  were  attached  to  copper  blocks  using  silver  epoxy  at 
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around  90°C.  The  test  results  are  given  in  the  next  section 
along  with  comparison  with  the  theory. 


2 . 3 Experimental  Results 


The  probe  test  on  the  cleaved  bars  reveals  that  the 
threshold  current  is  very  uniform  across  the  wafer.  The 
measured  P-I  characteristic  (power  vs.  current)  for  these 
shiny  contact  lasers  at  300K  is  plotted  in  Figure  2-7  for 
different  cavity  lengths  under  0 . 5jls  pulse  at  1kHz . The 
lasers  with  L=125|lm  produce  1 . 1W  power  at  current  of  2A,  and 
the  500flm  ones  produce  0 . 5W  at  2A.  These  power  levels  are 
among  the  highest  for  InGaAsP/InP  DH  lasers,  and  are  close  to 
those  obtained  from  state-of-the-art  separate  confinement 
heterostructure  (SCH)  single  quantum  well  (SQW)  lasers 
[Zor93]  . Figure  2-8  shows  the  inverse  of  the  slope 
efficiency  vs.  cavity  length.  The  slope  efficiency  T|s  is 

defined  as  the  relative  increase  in  the  optical  power  with 
the  input  current  in  the  P-I  curve  beyond  threshold. 


Tls  = 


Ap 


= Vc  Tld  = Vc  TJi 


In (1/R) 


Al  ' ” '*  atL  + In  (1/R) 

where  Vc  is  the  characteristic  voltage 
corresponding  to  the  bandgap)  . The  measured  T|s 
is  curve  fitted  as 


(2 .11) 

(the  voltage 
in  Figure  2-8 


T)s  1 = 1.278  + 0.00245L  (2.12) 
where  T|s  is  in  W/A  and  L in  |lm.  Comparing  eqs . (2.12)  with 
(2.11),  one  gets  at=22cm_1  and  T|i=82%, 


very  reasonable  values 
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in  good  agreement  with  the  theoretical  assumption  in  the 
previous  section.  Such  a low  value  for  at  indicates  that  the 

shiny  contact  introduces  negligible  optical  loss,  consistent 
with  the  theoretical  calculations  in  Section  2.1.  The 
threshold  current  is  plotted  in  Figure  2-9  as  a function  of 
cavity  length.  It  is  seen  that  the  experimental  data  is 
fitted  as  the  following  expression 

Ith  = 153.5  + 1.184L  (2.13) 
where  Ith  is  in  mA  and  L in  Jim.  Compare  this  equation  with 
the  theory. 


, T i r at  3tr,  , wdln(l/R) 

Ith  = Lwd  [ — + ] + - ' 


(2.14) 


Usp  P r TJgp  TJgpPr 

one  obtains  r|sp(3  = 0 . 0275cm-|fm/A  and  jtr/Tlsp  = 4 4 4 6 A/ cm2-Jlm  . 
Assume  the  calculated  value  for  jtr  is  3000A/cm2-^m,  then 
ilsp=57%  and  hence  Tir=82%.  The  Auger  coefficient  and  the  gain 
slope  should  be  C=1 . 6xl0_29cm6/s  and  P=0 . 0334cm-flm/A . One  can 
see  that  both  C and  (3  are  overestimated  in  the  theory  by 

about  88%  and  135%  respectively.  Both  the  theoretical  and 
experimental  results  for  C and  P are  in  the  same  order  of 

magnitude  and  are  reasonably  close,  considering  the 
uncertainties  in  the  energy  band  structure  for  InGaAsP 
[Kat  92 ] . 


The  P-I  curve  for  the  L=250|lm  laser  at  different 
temperatures  is  plotted  in  Figure  2-10,  and  the  temperature 
dependence  of  the  threshold  current  and  slope  efficiency 
given  in  Figures  2-11  and  2-12  respectively.  It  is  clear 
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that  the  threshold  current  and  slope  efficiency  can  be 
approximated  as 

Ith  = 300  exp (T/59)  (2.15) 

71s  = 0.641  - 0.00276T  (2.16) 
where  Ith?  Us  and  T are  in  mA,  W/A  and  °C  respectively.  The 
measured  characteristic  temperature  To=59°C  and  the 
temperature  dependence  of  Tls  agree  excellently  with  previous 
reports  [Agr86] . 

In  the  next  section  we  discuss  the  expected  improvements 
in  the  CW  performance  for  such  shiny  contact  lasers  relative 
to  conventional  DH  lasers  using  parameter  values  extracted 
from  the  measurements  presented  in  this  section. 

2 . 4 Further  Discussion 

In  Section  2.3  it  has  been  shown  that  the  pulsed 
performance  of  the  1 . 3p.m  InGaAsP/InP  shiny  contact  laser  is 
consistent  with  the  theoretical  predictions.  In  this 
section,  the  CW  performance  for  such  lasers  is  investigated 
at  room  temperature,  in  which  case  the  thermal  impedance  is  a 
key  parameter. 

According  to  the  analysis  presented  in  Section  2.1,  the 
optimal  p-clad  thickness  for  a 1 . 3|lm  InGaAsP/InP  shiny 
contact  laser  structure  should  be  t=0.5|lm.  In  applications 
where  single  mode  operation  is  preferred,  the  cavity  length 
and  the  window  width  should  not  be  too  large . Therefore  in 
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this  section  a shiny  contact  structure  with  the  following 
dimension  parameters  is  assumed:  t=0.5p.m,  L=250|J.m  and  w=4|lm. 

First  of  all,  the  output  power  vs.  current  for  a diode 
laser  can  be  expressed  as 

po  = Tls(I-Ith)  (2.17) 
where  the  threshold  current  I^h  and  the  slope  efficiency  TJS 
are  assumed  to  have  the  following  forms,  according  to  our 
measurements, 

1th  = 0.025  exp [ (Ta-Ts)  / 5 9 ] (in  A)  (2.18) 
232-T 

T[s  = 0-48  ^232-Tg^  (in  W/A)  (2.19) 
where  Ta  and  Ts  are  the  temperatures  in  °C  of  the  active  layer 
and  the  heat  sink  respectively  (TS=27°C)  . A lower  injection 
efficiency  (70%)  has  been  used  for  such  a narrow  stripe 
structure  to  account  for  the  lateral  current  flow.  The  I-V 
(current-voltage)  relation  of  a diode  laser  is 

V = Vc  + I rd  (2.20) 
where  VC=1V  and  rd=10£2  are  assumed.  Note  that  not  all  of  the 
input  electrical  power  is  emitted  as  light  but  a large 
fraction  of  it  is  converted  into  heat;  hence,  one  may  write 

Q=  VI  - P0  = (Ta-Tg)/Rth  (2.21) 
where  Q is  the  thermal  power  (heat  flux)  and  Rth  the  thermal 
impedance.  If  Rth  is  known,  eqs . (2  . 17 ) - ( 2 . 2 1 ) can  be  solved 
self  consistently  to  obtain  the  actual  temperature  in  the 
active  layer  Ta  at  a given  current,  from  which  the  output 
power  P0  is  then  calculated. 
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The  thermal  impedance  can  be  estimated  by  assuming  that 
most  of  the  heat  is  generated  in  the  center  of  the  active 
layer  and  flows  through  the  epitaxial  layers  to  the  heat  sink 
one  dimensionally  [Kre77] . In  this  case,  the  thermal 
impedance  of  each  layer  can  be  written  as  pth  d/A  where  pth  is 

the  thermal  resistivity,  d the  length  along  the  heat  flow  and 
A the  area  of  the  cross  section.  The  following  values  are 
used  in  the  calculations  [Kat92]:  pth=l  • 4 7 cm-K/W  for  InP, 
Pth  = 2 0 . 9 cm-K/W  for  In.53Ga.47As  and  pth= 2 2 . 5 cm-K/W  for 
In . 73Ga . 27As  . 6P  . 4 • The  copper  block  accounts  for  about  18K/W 
approximately  [Kre77].  The  experimentally  measured  data  for 
Rthf  however,  generally  agrees  with  those  obtained  using  a two 
dimensional  analysis  [Yan81],  which  shows  that  Rth  should  be 
about  one  half  of  that  calculated  one  dimensionally  for  a 
narrow  stripe  laser.  Therefore  we  estimate  Rth  to  be  40K/W 
for  the  shiny  contact  laser  (h=300A  and  t=0 . 5(lm)  and  100K/W 
for  the  conventional  laser  (h=0.5pm  and  t=l . 5pm) . 

The  predicted  CW  power  vs . current  is  plotted  in  Figure 
2-13.  Obviously,  both  the  linearity  and  highest  output  power 
are  shown  to  be  greatly  improved  using  the  shiny  contact 
laser  structure.  This  figure  confirms  that  in  conventional 
InGaAsP/InP  DH  laser  structures,  the  high  power  or  high 
temperature  performance  is  limited  by  the  heating  effect  due 
to  the  large  thermal  impedance.  For  example,  at  I=100mA,  the 
temperature  in  the  active  layer  is  43°C  for  the  conventional 
laser  while  it  is  33°C  for  the  shiny  contact  laser-- 
improvement  of  10°C  can  be  obtained.  At  higher  current 
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levels,  the  improvement  is  even  bigger:  Ta  is  81°C  vs.  47°C  at 
200mA  and  143°C  vs.  70°C  at  300mA. 

In  conclusion,  we  have  demonstrated  the  first  operation 
of  shiny  contact  lasers  in  the  1 . 3|lm  InGaAsP/InP  system,  and 
achieved  power  levels  greater  than  1W  pulsed  at  current  of  2A 
at  room  temperature  - comparable  to  those  obtained  from 
state-of-the-art  single  quantum  well  lasers.  The 
experimental  data  has  been  successfully  modeled  using  the 
theory,  and  important  parameters  extracted  from  the 
measurements.  The  improvements  in  the  CW  performance  for 
such  shiny  contact  lasers  relative  to  conventional  devices 
are  predicted  using  these  measured  parameters. 
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Figure  2-1:  The  InGaAsP/InP  DH  laser  structure 
used  in  the  theoretical  analysis 


optical  mode  loss  am  (cm-1) 
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Figure  2-2:  The  extra  optical  mode  loss 
for  the  structure  shown  in  Figure  2-1 


optical  mode  confinement  factor  r (%) 
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p-clad  thickness  t (jim) 


Figure  2-3:  The  optical  mode  confinement  factor 
for  the  structure  shown  in  Figure  2-1 


threshold  current  density  J,h  (A/cm2) 
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p-clad  thickness  t (j^m) 


Figure  2-4:  The  threshold  current  density 
for  the  structure  shown  in  Figure  2-1 
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Figure  2-5:  The  differential  quantum  efficiency 
for  the  structure  shown  in  Figure  2-1 


26 


w=100|lm 
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Figure  2-6:  A schematic  view  of  the  shiny  contact  laser 
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Figure  2-7:  The  measured  total  power  vs.  current  for  the 
InGaAsP  shiny  contact  lasers  fabricated  at  UF 


inverse  slope  efficiency  r|s_1  (A/W) 
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Figure  2-8:  The  inverse  slope  efficiency  vs.  cavity  length 
for  the  InGaAsP  shiny  contact  lasers  fabricated  at  UF 
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Figure  2-9:  The  threshold  current  vs.  cavity  length 
for  the  InGaAsP  shiny  contact  lasers  fabricated  at  UF 
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Figure  2-10:  The  measured  total  power  vs.  current 
for  the  L=250|im  laser  at  different  temperatures 


threshold  current  lth  (mA) 
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temperature  T (°C) 


Figure  2-11:  The  threshold  current  as  a function  of  temperature 
for  the  L=250|im  shiny  contact  laser 
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Figure  2-12:  The  slope  efficiency  as  a function  of  temperature 
for  the  L=250jlm  shiny  contact  laser 
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current  (mA) 


Figure  2-13:  The  calculated  CW  power  vs.  current  plot 
for  the  optimized  InGaAsP  shiny  contact  laser 


CHAPTER  III 

FLOQUET-BLOCH  ANALYSIS  OF  THE  SHINY  CONTACT 
DISTRIBUTED  FEEDBACK  DIODE  LASER 


In  the  previous  chapter,  the  operation  of  1 . 3}lm 
InGaAsP/InP  shiny  contact  lasers  has  been  demonstrated  for 


the  first 

time,  and 

several 

advantages 

of  such 

devices 

described . 

In  this 

chapter, 

we  discuss 

a shiny 

contact 

distributed  feedback  (DFB)  laser  structure  in  which  the 
diffraction  grating  is  a corrugated  metal/semiconductor 
interface  as  shown  in  Figure  3-1.  Heterostructure  diode 
lasers  with  such  integrated,  metalized  diffraction  gratings 
were  first  demonstrated  about  eighteen  years  ago  [Zor75, 
Wal77].  The  two  main  reasons  for  using  metalized  gratings 
instead  of  dielectric  gratings  are  thermal  impedance 
minimization  and  elimination  of  epitaxial  regrowth  after 
grating  fabrication.  The  former  has  been  explained  in 
Chapter  II,  and  the  latter  could  not  only  reduce  the 
fabrication  cost  but  also  improve  the  reliability.  In 
conventional  DFB  lasers,  the  diffraction  grating  is  only 
about  0 . l(i.m  from  the  active  layer  and  one  has  to  perform  a 
high  temperature  epitaxial  regrowth  process  after  grating 
fabrication.  The  high  temperature  regrowth  on  a corrugated 
interface  can  possibly  diffuse  the  surface  defects  into  the 
active  layer  and  hence  introduces  a reliability  problem.  In 
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the  shiny  contact  DFB  laser  configuration,  the  diffraction 
grating  is  placed  about  an  optical  wavelength  (dimension 
t=0 . 3-0 . 5|lm  in  Figure  3-1)  from  the  active  layer  and  the 
ohmic  contact  is  made  directly  on  top  of  the  grating,  thus 
the  fabrication  cost  is  reduced  and  the  reliability  issue  is 
avoided  since  no  epitaxial  regrowth  is  required. 

Other  advantages  of  using  such  metalized  gratings 
include  better  single  frequency  yield  for  edge  emitting  DFB 
lasers  [Sha89]  and  higher  CW  power  for  grating  surface 
emitting  (GSE)  lasers  [Mac87]  . The  former  is  due  to  the 
combined  index/loss  DFB  coupling  mechanism  built  in  the 
metalized  grating,  which  is  absent  in  conventional  type  DFB 
lasers  using  dielectric  gratings.  The  latter  is  because 
there  is  only  one  radiation  beam  reflecting  from  the 
metalized  grating,  in  contrast  with  two  radiation  beams  (one 
upward  and  one  downward)  from  conventional  GSE  lasers  with 
dielectric  gratings,  and  one  of  the  two  beams  is  not 
utilized . 

The  main  problem  with  metalized  gratings  is  that  the 
corrugated  metal/semiconductor  interface  must  have  low 
electrical  resistance  and  low  optical  mode  loss--shiny 
contact  grating.  Our  experiment  shows  that  this  problem  can 
be  solved  [Luo92]  and  recent  reports  indicate  that  high 
performance,  reliable  devices  have  been  fabricated  [Sha88, 
Mot 89,  Ras91 ] . 

In  modeling  laser  action  in  DFB  devices,  one  normally 
utilizes  coupled  mode  theory  [Kog72]  or  transfer  matrix 
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approach  [Han92].  These  standard  techniques  do  not  work  well 
in  the  shiny  contact  DFB  structure  due  to  the  large  index 
difference  at  such  metalized  gratings  (see  Chapter  IV) . As  a 
result,  a Floquet-Bloch  formalism  [Nol90]  is  developed  in 
this  chapter  to  analyze  the  wave  propagation  in  the  shiny 
contact  DFB  laser.  The  full  Floquet-Bloch  expansion  is 
truncated  by  including  only  the  minimum  number  of  terms  and 
applied  to  several  practical  configurations.  Numerical 
results  are  presented  in  the  next  chapter  for  comparison  with 
those  obtained  using  the  ray  optics  and  standard  approaches. 

3 . 1 Full  Floquet-Bloch  Formalism 

The  configuration  used  in  the  analysis  is  that  of  a 
perfect  metal  clad  sinusoidally  corrugated  waveguide  for  TE 
modes,  i.e.,  the  electric  field  is  y-oriented  in  Figure  3-1. 
Assume  the  electric  field  mode  in  each  layer  is  expanded 
using  the  following  plane  wave  series  based  on  Floquet-Bloch 
theorem: 

oo 

Ed(x,z)  = Amexp [-qm(x-t-d) ] exp (-iPmz) 

m=-oo 

oo 

Ea(x,z)  = ^ { Bmcos  [Gm(x-t)  ] +Cmsin  [Om(x-t)  ] } exp  ( -i(3mz ) 

m=-oo 

oo 

Eb(x,z)  = X { Dmexp  (pmx) +Emexp  (-pmx)  } exp  (-i^z) 

m=-oo 

Ec(x,z)  = 0 in  the  perfect  metal  (3.1) 


where 
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Pm  = Po  + mK  = Pq  + m(27C/A) 


Pm  = VPm2-nb2k2 


qm  = Vpm2-nd2k2 


k = 2%/X 


X:  wavelength  in  vacuum 


Rayleigh  was  the  first  to  use  a similar  plane  wave  series  in 
analyzing  classical  diffraction  problems  and  he  assumed  that 
such  a plane  wave  series  can  uniformly  converge  to  the  unique 
solution  when  the  number  of  terms  included  goes  to  infinity, 
and  this  assumption  is  referred  as  Rayleigh  assumption 
[Pet80]  . Despite  the  fact  that  Rayleigh  assumption  is  not 
always  true  for  an  arbitrary  grating  shape,  the  results 
obtained  using  this  expansion  are  generally  fast  and  accurate 
for  sinusoidal  groove  gratings  [Pet80]  . Consequently  this 
method  is  adopted  to  investigate  the  wave  propagation  in  the 
shiny  contact  DFB  structure. 

At  x=t,  the  electric  field  and  its  derivative  should  be 
continuous,  and  this  leads  to 


oo  oo 


X Bmexp(-iPmZ)  = X i DmexP <Pmt ) +Emexp (-pmt ) } exp (-iPmz) 


Bm  = Dmexp  (pmt ) +Emexp  ( -pmt ) Vm 


(3.2) 


oo  OO 


Similarly  at  x=t+d. 
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oo  oo 

X Amexp(-i(3mz)  = £ { Bmcos  (amd) +Cmsin  (omd)  } exp  (-ipmz) 

m=-oo  m=-°° 

=>  Am  = Bmcos  (amd) +Cmsin  (amd)  Vm  (3.4) 

OO  OO 

X cImAme xp  ( - i PmZ ) = Gm  { Bmsin  (Omd)  -Cmcos  (Omd)  } exp  ( -iJ3mz ) 

m=-oo  m=-°o 

=>  qmAm  = Om{  Bmsin  (amd) -Cmcos  (Omd)  } Vm  (3.5) 

At  the  grating  surface  x=f (z) =acos (Kz) , should  vanish 

since  a perfect  metal  has  been  assumed, 

OO 

X { Dmexp [pmf (z) ] +Emexp [-pmf (z) ] } exp  (-i^mZ)  =0  (3.6) 

m=-oo 

Recall  the  following  Fourier  series, 

OO 

exp [±pmacos (Kz) ] = ^ In (±pma)  exp (inKz) 

n=-°o 

where  In(u)  is  the  modified  (or  hyperbolic)  Bessel  function 
of  the  first  kind  (of  order  n) . The  following  equations  are 
listed  for  reference. 

In(-u)=(-l)nIn(u) 

I-n  (u)  =In  (u) 

Then  eq.(3.6)  becomes 

00  oo 

X exp(-iPmz)  { Dm  £ In (pma) exp (inKz) 
m=— oo  n=-o° 

OO 

+ Em  X (-1)  nIn  (pma)  exp  (inKz)  } = 0 

n=-oo 
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X exp(-if3oz)  { Dm  X In  (pma)  exp  [-i  (m-n)  Kz] 

m=-oo  n=-oo 

oo 

+ Em  X (-1)  nIn  (pma)  exp  [-i  (m-n)  Kz]  } = 0 

n=-oo 

Since  the  coefficient  of  exp(-inKz)  for  each  n should  be 
zero,  one  then  obtains 

CO 

^ { Dmim-n  (Pma ) + Em(-1)  lm-n(Pma)  } = 0 Vn  (3.7) 

m=-oo 

Equations  (3.2)- (3.5)  and  (3.7)  constitute  an  infinite  linear 
homogeneous  system.  The  number  of  equations  is  equal  to  the 
number  of  variables,  and  therefore  a nonzero  solution 
requires  that  the  determinant  vanish.  Thus  a characteristic 
eigenvalue  equation  can  be  obtained  and  solved  for  Po  at  each 
A,--dispersion  relation. 

In  practice  only  finite  terms  can  be  considered  in  this 
linear  system.  Suppose  one  keeps  N orders  (there  are  N 
values  for  the  order  index  m) , then  one  gets  a 5NX5N  matrix! 

Due  to  the  huge  amount  of  numerical  computation  involved,  it 
is  very  difficult  to  present  a general  treatment. 
Consequently,  several  special  cases  of  practical  interest  are 
studied  based  on  a simplified,  truncated  expansion  in  the 
next  section.  The  numerical  results  obtained  using  the 
truncated  Floquet-Bloch  formalism  are  given  in  Chapter  IV 
along  with  those  from  the  ray  optics  and  other  standard 
techniques  for  comparison. 
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3 . 2 Truncated  Floquet-Bloch  Formalism 

First  of  all  we  look  at  the  zeroth  order  case  where  only 
m=0  is  considered.  This  case  represents  an  unperturbed, 
planar  four-layer  waveguide  with  a perfect  metal  clad,  i.e., 
a=0 . The  above  linear  homogeneous  system  leads  to 
Bo  - Doexp(pot)  - Eoexp(-pot)  = 0 
CT0C0  - PoDoexp(p0t)  + p0E0exp  (-p0t)  = 0 
Aq  - Bqcos (Ood)  - Cosin (Ood)  = 0 
qoA0  - aoBosin(Ood)  + OqCoCOS (O0d)  = 0 
Do  + Eo  = 0 

where  Io(0)=l  and  In(0)=0  for  n^O  have  been  used.  Requiring  a 
nontrivial  solution,  this  5x5  system  can  be  manipulated  to 

obtain  the  following  transcendental  equation: 

pocoth(pot)  = Ootan [Ood-tan-1 (q0/ao)  ] (3.8) 

Due  to  the  mirror  principle,  this  equation  is  the  same  as  for 
the  odd  mode  of  a symmetric  5-layer  waveguide  whose 
refractive  index  is  nd/na/nb/na/nd  and  thickness  oo/d/2t/d/°° . 
The  dispersion  curve  for  such  a planar  waveguide  is  shown  as 
the  solid  line  in  Figure  3-2 . 

Next  let  us  consider  a first  order  grating  with  a 
shallow  groove,  i.e.,  K=2po  and  ka«l.  In  this  case,  there 

are  only  two  traveling  terms  that  carry  electromagnetic 
power:  the  fundamental  forward  wave  (m=0,  propagation 

constant  = Po)  and  the  backward  diffracted  wave  (m=-l, 
propagation  const ant=-po ) • All  higher  order  waves  are 
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evanescent  along  x and  fast  oscillating  along  z,  and  are 
neglected  in  the  following  calculations.  This  two  wave 
approximation  is  valid  for  shallow  groove  gratings  normally 
utilized  in  most  practical  DFB  laser  configurations.  From 
eqs . (3.2)-(3.5),  one  has 

Bq  - Doexp(pot)  - Eoexp(-pot)  = 0 

B-i  - D-iexp(p-it)  - E-iexp (-p_it)  = 0 

<*oCo  “ P0D0exP  (P0t ) + PoE0exp  (-p0t)  = 0 

a-iC-i  - p-iD_!exp (p-it)  + p-iE-xexp (-p_!t ) = 0 

Aq  - Bqcos  (dod)  - Cosin  (God)  = 0 

A_i  - B_icos (O-id)  - C-isin(G-id)  = 0 

qoAo  - GoBosin(God)  + GoCqCOS (God)  = 0 


q_lA_i  - G_iB_isin  (G_id)  + G_iC_icos (G_id)  = 0 
In  eq.  (3.7),  one  has  to  remember  that  the  Bessel  function  of 
small  argument  is  fast  decaying  with  increasing  order.  Hence 
one  needs  to  pick  up  the  most  significant  two  equations  from 
eq.  (3.7) . In  the  present  case  for  first  order  gratings,  one 
should  choose  n=0  and  n=-l: 


DoIo(Poa)  + Eolo(poa)  + D-iI-i(p-ia)  - E-iI-i (p-ia)  = 0 
Doli(poa)  - Eoli(poa)  + D-iIo(p-ia)  + E-iIo(p-ia)  = 0 
After  rearrangement,  this  10x10  matrix  leads  to  the  following 
eigenvalue  equation: 

Ip (ppa) • Ip (p-ia) 


II (p0a) • Ii (p_ia) 


(3.9) 


where 
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tanh  { tanh  1 { (— ) tan  [Oid+tan-1  (— ) ] } + pit  } 

Oi  li 

Since  we  deal  with  shallow  gratings  here,  the  Bessel 
functions  in  eq.  (3.9)  can  be  further  approximated  as 
polynomials : 

Io(u)  = 1 

for  u«l 

II (u)  = u/2 

The  dispersion  curve  for  such  a corrugated  waveguide  with  a 
first  order  grating  is  shown  as  the  thick  dashed  line  in 
Figure  3.2.  It  is  seen  that  without  the  corrugation 

(perturbation)  , the  CO — Po  dispersion  curve  is  continuous  and 

smooth.  With  the  perturbation  turned  on,  the  interaction 
between  the  waveguide  mode  and  the  grating  introduces  a stop 
band  Ago  in  the  CO— po  curve  at  the  Bragg  frequency  Jt/A.  Inside 
the  stop  band,  the  propagation  constant  Po  becomes  a complex 
number  and  the  forward  wave  is  attenuated  along  the 
corrugated  waveguide.  This  attenuation  is  due  to  the 
continuous  feedback  under  phase  matched  condition  rather  than 
absorption  or  loss,  i.e.,  the  power  of  the  forward  wave  is 
continuously  reflected  into  that  of  the  backward  wave.  The 
stop  band  in  Figure  3-2  is  very  similar  to  the  bandgap  in  the 
energy-wavevector  diagram  for  the  electron  state  in 
semiconductors  . 

The  stop  band  of  the  dispersion  curve  describes  the 
strength  of  the  interaction  between  the  grating  and  the 
optical  mode,  and  its  magnitude  is  related  to  the  coupling 
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coefficient  K (introduced  in  the  next  chapter)  according  to 

coupled  mode  theory  [Yar73]  : 

2 K c 

A co  = 

nef  f 

AX 

K = 71  neff  2 (3.10) 

where  c is  the  velocity  of  light  in  vacuum  and  neff  the 
effective  waveguide  index  for  the  transverse  mode  of 
interest.  Equation  (3.10)  is  used  in  Chapter  IV  to  verify 
the  accuracy  of  the  ray  optics  expression  for  K. 

For  a second  order  grating  with  a shallow  groove,  K=(3q 

and  ka«l,  there  are  now  three  propagating  waves:  the  forward 
going  wave  (m=0,  propagation  constant=(3o)  t the  upward  surface 

emitting  wave  (m=-l,  propagation  constant=0)  and  the  backward 
diffracted  wave  (m=-2,  propagation  constant  = -(3o ) . As  a 

result,  one  can  obtain  a 15x15  matrix  from  the  homogeneous 

system,  eqs . (3.2)-(3.5)  and  (3.7)  by  neglecting  all 

evanescent  waves.  Following  a similar  procedure,  the  15x15 

matrix  leads  to  the  transcendental  characteristic  equation: 

I0  (p-ia)  [Io  (poa)  Io  (P-2a)  -I2  (Poa)  I2  (p-2a)  ] 

= Ii  (p_ia)  ¥-2  { Ii  (poa)  [Io  (P-2&)  -I2  (P-2a)  ] ^0 

+ Ii (p-2a)  [Io (poa) -I2 (Poa) ] ^-2  } (3.11) 

Again  eq.  (3.11)  is  valid  only  for  shallow  groove  gratings 
because  all  the  higher  order  evanescent  waves  are  neglected 
in  the  derivation.  The  Bessel  functions  can  be  approximated 
as  polynomials  for  small  argument : 
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I0(u)  = 1 + ~f~ 

II  (u)  = ~~  for  u«l 

I2(u)  = ~- 


In  the  next  chapter,  the  ray  optics  formalism  is 
introduced  to  derive  closed-form  expressions  for  the  coupling 
coefficient  K in  conventional  type  and  shiny  contact  DFB 

laser  configurations.  Numerical  results  from  the  Floquet- 
Bloch,  ray  optics  and  other  standard  approaches  are  compared 
in  various  DFB  laser  configurations.  The  advantages  and 
disadvantages  of  each  approach  are  discussed. 


2.  active 


x=d+t 

x=t 


n. 


t 


3.  buffer  nb  t 


4.  perfect  metal  nc=-i°o 


grating  function:  x = f(z)  = a-cos(Kz)  = a • cos(27tz/A) 


Figure  3-1:  The  shiny  contact  DFB  laser  structure 
used  in  the  analysis 
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propagation  constant  Re(p0) 


Figure  3-2:  The  co — Pq  dispersion  relation 
for  the  structure  shown  in  Figure  3-1 


CHAPTER  IV 
RAY  OPTICS  APPROACH 

TO  THE  DISTRIBUTED  FEEDBACK  DIODE  LASER 


Distributed  feedback  ( DFB ) diode  lasers  are  key 
components  in  modern  optical  communication  systems.  In 
modeling  laser  action  in  these  devices,  one  normally  utilizes 
coupled  mode  theory  [Kog72]  . In  this  theory,  the  key  laser 
parameter  is  the  backward  coupling  coefficient  K since  its 

magnitude  determines  such  parameters  as  frequency  selectivity 
and  threshold  mode  gain.  The  standard  technique  to  determine 
K in  a DFB  laser  is  to  evaluate  the  coupled  mode  integral 

which  describes  the  overlap  of  the  unperturbed  mode  with  the 
grating  perturbation  [Str75]: 

k 2 j*  * 2 2 

K = — An  (x)  E^  (x)  dx  (4.1) 

2p  J 

where  k is  the  wavevector  in  vacuum,  (3  propagation  constant 

o 

of  the  waveguide  mode,  An  (x)  index  perturbation  and  E (x) 
normalized  electric  field  mode. 

This  technique  works  well  for  weak  perturbations  but 
shows  an  inherent  problem  as  to  how  to  choose  the  unperturbed 
waveguide  geometry  when  the  perturbation  is  strong  [Str77, 
Cor88] . For  a corrugated  semiconduct or/per feet  metal 
interface,  the  coupled-mode  overlap  integral  cannot  be 
evaluated  at  all  since  the  index  perturbation  is  infinite 
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while  the  electric  field  is  zero  and  hence  the  product  is 
undefined.  For  a realistic  semiconduct or/metal  or 
semiconductor/air  grating  where  the  index  difference  is 
large,  the  standard  approach  gives  results  strongly  dependent 
on  the  choice  of  the  unperturbed  waveguide  [Sha89] . 

Another  widely  used  approach  in  modeling  DFB  lasers  is 
the  effective  index  transfer  matrix  method  [Han92].  In  this 
approach  one  transforms  the  corrugated  waveguide  into  a stack 
of  homogeneous  layers  with  local  effective  indices.  After 
the  transformation,  a standard  matrix  formalism  from 
classical  optics  is  used  to  model  the  spectral  properties. 
For  comparison  with  coupled-mode  theory,  a definition  for  the 
coupling  coefficient  is  given  [Har84]  . In  contrast  to 
coupled  mode  theory,  the  transfer  matrix  method  does  not 
involve  an  unperturbed  structure  and  hence  avoids  the 
uncertainty  problem.  However,  there  are  difficulties  in 
considering  transverse  modes  and  ambiguities  in  such 
transformations  for  highly  asymmetrical  corrugated 
waveguides.  For  example,  in  a simple  three  layer  DFB 
waveguide  with  a metal  clad  on  one  side  and  a dielectric  clad 
on  the  other  side,  the  transfer  matrix  method  gives  exactly 
the  same  results  for  whether  the  grating  is  placed  on  the 
metal  side  or  the  dielectric  side,  i.e.,  such  transformations 
do  not  distinguish  the  two  different  cases.  Therefore  one 
must  be  careful  when  representing  a corrugated  waveguide  by  a 
stack  of  homogeneous  layers — the  standard  transformation 
should  be  adjusted  according  to  some  known  reference. 
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In  this  chapter,  a ray  optics  approach  to  the  first 
order  DFB  coupling  coefficient  is  presented,  and  closed-form 
expressions  for  K are  derived  in  DFB  configurations  with 

metalized  and  dielectric  gratings.  For  standard  type  DFB 
lasers  with  dielectric  gratings  where  the  index  difference  is 
small,  the  K values  from  the  coupled-mode  overlap  integral, 

the  transfer  matrix  and  the  ray  optics  methods  agree  very 
well.  For  a distributed  Bragg  reflector  (DBR)  type  structure 
with  a corrugated  semiconductor/air  interface  where  the  index 
difference  is  larger,  the  K values  from  the  ray  optics  method 

agree  with  those  from  the  transfer  matrix  method,  but  the 
coupled-mode  overlap  integral  gives  inconsistently  smaller 
results  depending  on  the  choice  of  the  unperturbed  waveguide. 
For  shiny  contact  DFB  lasers  with  metalized  gratings  where 
the  index  difference  is  very  large,  the  standard  techniques 
do  not  work  well  and  the  accuracy  of  the  ray  optics 
expression  is  verified  by  comparing  the  results  with  those 
obtained  using  the  truncated  Floquet-Bloch  approach  explained 
in  Chapter  III.  The  ray  optics  formalism  is  further  extended 
in  Chapter  V to  deep  groove  configurations  where  the  transfer 
matrix  fails. 


4 . 1 General  Formalism 


The  DFB  coupling  coefficient  K is  a measure  of  the 
backward  coupling  strength  per  unit  length,  so  from  a ray 


50 


optics  point  of  view,  K can  be  expressed  by  the  following 
equation  [Luo90,  Luo92] : 


where  B is  the  bounce  rate  for  the  transverse  mode  of 
interest  in  the  planarized  version  of  the  waveguide  and  T]  the 

diffraction  efficiency  of  the  grating  at  mode  bounce  angle 
0eff  (see  Figure  4-1) . The  bounce  rate  is  given  by 


where  na  is  the  refractive  index  of  the  active  layer, 
nef f=nasin0eff  the  effective  refractive  index  for  the  mode, 
and  wef f the  effective  waveguide  thickness.  The  effective 
waveguide  thickness  weff  is  a measure  of  the  decay  of  the 
electric  field  mode  outside  the  active  layer  and  is  related 
to  the  optical  confinement  factor.  The  expression  for  weff 
can  be  derived  directly  from  the  definition  [Tam90]  for  each 
waveguide  structure,  and  is  given  in  the  following  sections. 

A closed-form  expression  for  the  diffraction  efficiency 
can  be  derived  for  shallow  grooves  (ka«l)  in  two  ways.  The 
first  method  is  very  similar  to  the  truncated  Floquet-Bloch 
approach  as  presented  in  Section  3.2.  Assume  a TE  polarized 
plane  wave  is  incident  at  the  na/nb  interface  at  angle  6eff  as 
shown  in  Figure  4-1. 

Ei  = Aq  exp(-iGx  - i(3z) 


(4.2) 


(4.3) 


2 neff  weff 
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where  P=kneff  is  the  propagation  constant  and  a=k^na2  - neff2. 
The  reflected  waves  resulting  from  the  zero  and  first  order 
diffraction  can  be  written  as 

Er*0'  = ao  exp  (iax  - iPz ) 

Er(-1)  = a_i  exp  (iax  + i(3z) 

where  the  grating  periodicity  has  been  assumed  to  be  equal  to 
X,/2neff  since  we  deal  with  first  order  gratings.  Neglecting 

higher  order  waves  (see  Section  3.2),  the  electric  field  in 
the  incident  layer  na  is  then 

Ea  = Ei  + Er(0)  + Er(_1) 

= A0  exp  (-iax  - iPz) 

+ ao  exp (iax  - iPz)  + a_i  exp ( iax  + iPz ) (4.4) 

Similarly  for  layers  nb  and  nc,  the  total  electric  fields  can 
be  written  as 


Eb  = Bo  exp  (-px  - iPz)  + B_i  exp  (-px  + ipz) 


+ bo  exp  (px  - iPz)  + b_i  exp  (px  + iPz) 
Ec  = Co  exp  (-qx  - ipz)  + C_i  exp  (-qx  + iPz) 


(4.5) 

(4.6) 


where  p=k-\j  neff2  - nb2  and  q=k^  neff2  - nc2  . In  eqs  . (4.  4) -(4. 6), 

there  are  8 unknown  amplitudes  relative  to  the  incident 
amplitude  Ao,  and  therefore  one  needs  8 independent  equations 
(8x8  linear  inhomogeneous  system)  in  order  to  solve  the  first 


order  diffraction  efficiency: 

a-i  i 2 


T1  = 


A0 


2 


(4.7) 
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Matching  the  boundary  conditions  at  x=0  and  x=t+f(z)  results 
in  4 equations,  each  of  which  can  be  separated  into  two 
equations  by  setting  the  coefficients  of  exp(±i(3z)  to  zero. 
Therefore  a closed-form  expression  for  T|  for  ka<<l  can  be 
obtained . 

The  second  method  is  based  on  a bouncing  ray  picture  and 
is  physically  appealing.  Assume  a TE  polarized  plane  wave  is 
incident  at  the  na/nb  interface  at  angle  0eff-  Each  time  the 

ray  hits  the  planar  na/nb  interface,  a transmitted  and 
reflected  rays  are  generated  (denoted  by  tab,  tba  and  rab,  rba 

respectively)  . Each  time  the  ray  is  incident  at  the 
corrugated  nb/nc  interface,  there  are  four  rays  produced:  tb°* 
and  rb°’  are  the  0th  order  transmitted  and  reflected  rays,  and 
t^1’  and  are  the  1st  order  transmitted  and  reflected 
rays.  For  shallow  gratings  of  ka«l,  both  tj^11  and  r^1*  are 

much  smaller  than  1,  and  hence  in  the  following  derivation 
any  product  of  two  or  more  1st  order  terms  are  neglected.  The 
contribution  to  T|  includes  these  terms: 


} 


This  leads  to  a series  for  T) 


OO 


^ab^ba  rbc 


n=0 


(4.8) 
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This  series  is  simply  a derivative  of  a geometric  series,  and 
therefore  the  sum  is  in  a simple  form: 


VtT  = 


t t r^e'2^ 
^ab^ba  rbc  e 

n+r  r <0>  e-2^)  2 
' 1+rab  rbc  e ' 


(4.9) 


Expressions  for  r^°’  and  r^1’  are  a classical  problem  and  can 


be  obtained  from  [Mar72]  for  shallow  sinusoidal  grooves. 


(0) 


-Erg 


(-it  . p-g 

rbc  = 


(4.10) 


fc>c  p+g  *r'“  p+g 

Other  transmission  and  reflection  coefficients  are  written 
down  here  for  reference. 


20 


'ab 


rab  = 


O-ip 

G+ip 


'ba 


-2  ip 
O-ip 


(4.11) 


= -r 


ba 


O-ip 

Combining  eqs  . (4.9),  (4.2)  and  (4.3),  one  obtains  the 

following  expression  for  K, 


K = 


Vna2  ' n 


eff 


t t r( 
nabrba  rbc  e 


2 neff  weff  n+r  r(0)P-2ptv 
u+rabrbce  ’ 


(4.12) 


Note  that  the  above  two  methods  are  general  for  any  type 
of  cladding  layer  as  long  as  the  proper  refractive  index  is 
used,  e.g.,  nc  is  complex  for  a metalized  or  a 
gain/absorption  grating.  It  can  be  shown  that  eq.  (4.7) 
agrees  exactly  with  eq.  (4.9)  since  they  are  both  derived 
based  on  the  same  assumption  ka<<l . This  is  similar  to  the 
fact  that  in  a Fabry-Perot  resonator,  the  wave  optics  and  the 
ray  optics  lead  to  the  same  results. 
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It  is  interesting  that  the  unperturbed  waveguide 
uncertainty  problem  is  avoided  using  the  ray  optics  technique 
since  the  key  parameter  T]  is  derived  directly  from  the 

corrugated  waveguide  structure.  The  ray  optics  technique 
works  well  for  strong  coupling  as  well  as  weak  coupling  DFB 
configurations,  and  it  also  provides  device  designers  with 
simple  expressions  for  K which  are  much  easier  to  implement 

than  the  standard  overlap  integral  programs  since  they  can  be 
evaluated  without  using  numerical  techniques  of  any  kind  (see 
Section  4.3)  . In  addition,  the  ray  optics  technique  is 
physically  intuitive  and  gives  one  a different  perspective  on 
how  DFB  lasers  actually  work.  In  the  following  sections  the 
ray  optics  technique  is  applied  to  several  conventional  and 
nonconvent ional , including  the  shiny  contact  DFB  laser 
configurations . 

4 . 2 DFB  Lasers  with  Metalized  Gratings 

In  this  section,  we  show  how  to  derive  an  accurate 
closed-form  expression  for  K for  the  shiny  contact  DFB 

waveguide  shown  in  Figure  3-1  despite  the  complexity  of  such 
structures.  The  approach  we  use  to  demonstrate  the  accuracy 
of  the  ray  optics  technique  is  to  compare  ray  optics  and 
Floquet-Bloch  derived  plots  of  K for  perfect  metal  grating 

(see  Chapter  III)  . We  then  show  that  the  ray  optics  plots 
change  very  little  when  the  perfect  metal  is  replaced  by  a 
shiny  metal  such  as  gold. 
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A perfect  metal  is  by  definition  a metal  whose 

conductivity  is  infinite,  i.e.,  its  complex  refraction  index 
nc=-i°°.  In  this  case  q=°°,  r^°)=-l  and  r^J^ipa,  and  '\[r\  and 

weff  are 


= 


( ka ) V na2-neff2  ( neff2-nb2) 


(na-neff  ) sinlr  (pt)  + (neff2-nb2)  coshMpt ) 


1 (p2+o2) tanh (pt)  - G2pt/cosh2 (pt ) i 

W^ff  d + o o o t 

P p + a2tanh2(pt)  s 


1 1 
= d + — tanh(pt)  + 

p ^ s 


r 22 

where  s=kyneff  - nd  . The  coupling  coefficient  is  then,  from 


eq.  (4.12), 


K = 


ka  (na2-neff2)  (neff2-nb2) 


(4.13) 


2 neffweff  [ (na2-neff2)  sinh2  (pt)  + (neff2-nb2)  cosh2  (pt)  ] 

This  equation  is  useful  since  it  depends  only  on  the 
wavelength,  the  groove  depth  and  the  parameters  of  the  planar 
waveguide  structure.  Therefore,  once  the  eigenvalue  equation 
of  the  planar  waveguide  is  solved,  the  magnitude  of  the 
coupling  coefficient  can  be  calculated  in  a straightforward 
manner.  Strictly  speaking,  this  equation  is  valid  only  for 
ka<<l  because  of  the  assumption  needed  in  deriving  the 
expression  for  the  diffraction  efficiency.  However,  as 
indicated  by  Rigrod  and  Marcuse  [Rig76],  expressions  of  this 
type  are  in  most  cases  at  least  qualitatively  correct  for  ka 
values  well  beyond  the  ka«l  limit  used  in  the  derivation. 

In  the  real  world,  perfect  metals  are  closely 
approximated  by  "shiny  metals,"  i.e.,  those  metals  with  a big 
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extinction  coefficient.  Gold  is  one  such  shiny  metal  at 
wavelength  of  0.85|Im  since  its  extinction  index  is  much 
greater  than  its  refractive  index  (nc=0 . 1 6-i5 . 3 ) . For  such 
realistic  metals,  one  can  simply  use  eq.  (4.13)  with  the 
proper  complex  refractive  index. 

As  a first  example,  Figure  4-2  shows  a plot  of  K vs . ka 

for  the  AlGaAs  corrugated  double  heterostructure . It  is 
observed  that  for  perfect  metal  gratings,  the  ray  optics  and 
the  Floquet-Bloch  curves  agree  very  well  at  small  values  of 
ka,  deviating  by  only  10%  at  ka=l . The  dashed  line  (for 
gold)  follows  the  solid  line  (for  perfect  metal)  closely  as 
expected.  In  Figure  4-3  we  plot  the  coupling  coefficient 
versus  the  thickness  of  the  buffer  layer  t.  Both  approaches 
are  well  consistent  over  the  whole  range  of  t . The  roll-over 
behavior  is  attributed  to  the  fact  that  as  t becomes  very 
small,  the  mode  is  pushed  far  away  from  the  grating  where  the 
interaction  takes  place,  causing  the  mode  overlap  to  decrease 
as  t is  reduced  further.  In  Figure  4-4  the  coupling 
coefficient  is  plotted  versus  the  thickness  of  the  active 
layer  d,  with  t as  a parameter.  Excellent  agreement  between 
the  two  approaches  is  observed  again. 

The  above  derivation  is  based  on  gratings  of  a 
sinusoidal  shape.  For  other  grating  tooth  shape,  the 
coupling  coefficient  should  be  multiplied  by  a geometric 
factor  G:  G=1  for  sinusoidal  gratings,  G=8/tc2  for  triangular 
gratings  and  G=4/7C  for  square  gratings.  This  factor  comes 
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from  the  Fourier  coefficients  in  the  derivation  of  the 
diffraction  efficiency. 


In  this  section,  we  apply  the  ray  optics  technique  to  a 
generic  5-layer  separate  confinement  heterostructure  (SCH) 
DFB  laser  (see  Figure  4-5  (a)  ) , and  show  how  to  utilize  the 
ray  optics  expression  to  compute  K in  a more  complex, 

multiquantum  well  (MQW)  DFB  configuration  (see  Figure  4- 
5 (b)  ) . Then  the  results  from  the  ray  optics  expression  are 
compared  with  published  values  of  K determined  by  the 

standard  overlap  integral  and  transfer  matrix  techniques. 

The  effective  waveguide  thickness  weff  for  the  5-layer 
SCH  structure  is  given  as 


4 . 3 DFB  Lasers  with  Dielectric  Gratings 


weff  = d + 


(G2-p2)XY  + g2p(p2+q2)  (t  + l/q) 
p (O2  X2  + p2  Y2) 


+ 


(G2-p 1 2 ) X 1 Y 1 + P2p 1 (p 1 2+q 1 2 ) (t'+l/q') 

p'(o2x,2  + p,2Y'2) 


d + t + t ' + 


1 1 


where 


X = pcos(pt)  +qsin(pt) 


X'  = p'cos(p't')  +q'sin(p't') 
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Y = psin(pt)  -qcos(pt) 
After  manipulation,  eq.  (4.9) 
aap2  (p2+q2 ) 


a2x2+p2Y2 


Y'  = p'sin(p't')  -q'cos(p,t') 
leads  to 


aop2 (p2+q2) 

O2 [pcos (pt) +qsin  (pt) ] 2 + p2 [psin (pt) -qcos  (pt)  ] 2 

where  G is  the  geometric  factor  depending  on  the  groove  shape 
of  the  grating:  G=1  for  a sinusoidal  grating,  G=8/Jl2  for  a 

triangular  grating  and  G=4/k  for  a square  grating.  The 
coupling  coefficient  K is  then 


ka (na2-nef f2 ) p2 (p2+q2) 

K = G ' „ 0 0 — (4.14) 

2 neffWeff (02X2+p2Y2) 

g ka  (na2-neff2)  p2  (p2+q2) 

2neffweff  { O2  [pcos  (pt)  +qsin  (pt)  ] 2 + p2  [psin  (pt)  -qcos  (pt)  ] 2 } 

As  mentioned  above,  eq.  (4.14)  is  valid  only  for  ka<<l 
because  of  the  approximations  made  in  deriving  T] . However, 

as  indicated  by  Rigrod  and  Marcuse  [Rig76],  expressions  of 
this  type  are  in  most  cases  at  least  qualitatively  correct 
for  ka  values  well  beyond  the  limit  used  in  the  derivation. 

To  use  eq.  (4.14)  to  determine  K in  an  MQW  DFB  laser 

structure  as  shown  in  Figure  4-5 (b) , one  needs  to  transform 

the  MQW  structure  into  an  equivalent  SCH  structure  by 

replacing  the  quantum  well/barrier  layers  by  a single 

homogeneous  layer  with  an  r.m.s.  refractive  index  na  [Sam90] , 

2 2 
2 _ Nqwd-gwngw  + Nbrdbrnbr 

3 Nqwdqw  + Nbrdbr 


(4.15) 
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where  Nqw,  Nfc>r,  dqw , dbr/  nqw  and  n^r  are  the  numbers, 
thicknesses  and  refractive  indices  of  the  quantum  well  and 
barrier  layers  respectively. 

The  formula  for  K,  eq.  (4.14),  is  easy  to  use  since  it 
depends  only  on  the  wavelength,  the  grating  amplitude  and 
shape,  and  parameters  of  the  planar  waveguide  geometry. 
Therefore,  once  the  effective  index  neff  of  the  planar 
waveguide  is  obtained,  the  magnitude  of  the  coupling 
coefficient  can  be  calculated  in  a straightforward  manner. 
Since  the  accuracy  for  nm  is  not  critical  in  the  ray  optics 
computation,  one  can  use  an  approximate  method  to  estimate  nm 
without  solving  the  transcendental  eigenvalue  equation.  For 
example,  in  most  laser  structures,  nc=nc  and  the  thickness  of 
the  inner  layers  is  small  compared  with  the  wavelength. 
Therefore,  one  can  transform  the  SCH  or  MQW  waveguide  into  a 
simple  3 layer  waveguide  using  a technique  similar  to  eq. 
(4.15),  i.e.,  all  the  inner  layers  are  regarded  as  a single 
homogeneous  layer  with  an  r.m.s.  refractive  index  ne  and 
thickness  de . In  this  case,  ne  and  de  are  given  by 


and  the  following  closed-form  expression  for  neff  for  such  a 
symmetric  3-layer  waveguide  can  be  used  [Che83] , 


2 = dna2  + tnh2  + t 1 n^2 
d + t + t ' 


(4.16) 


de  = d + t + t ' 


(4.17) 


nef  f 


2 


2 


2 


2 


(4.18) 
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There  are  also  formulas  for  asymmetric  3-layer  waveguides 
(nc?tnc)  in  [Bot83]  and  the  references.  With  the  help  of 

these  formulas,  the  ray  optics  technique  provides  the  means 
for  quick  estimates  of  K in  relatively  complex  DFB 

configurations  without  using  numerical  techniques  of  any 
kind . 

The  above  derivation  is  written  for  neff<nb  and  neff<n^, 
i.e.,  the  transverse  components  of  the  two  counter  running 
waves  in  the  guiding  layers  are  traveling  waves.  In  case 
they  are  evanescent,  i.e.,  neff>nb  and  neff>n^,  then  p and  p' 
should  be  replaced  by  -iki/ neff2-nb2  and  -ik^  neff2-n&2 
respectively . 

As  a first  example,  we  plot  K versus  ka  in  Figure  4-6 
for  a typical  1.55|lm  symmetric  SCH  DFB  laser  [Sak82].  It  is 
clear  that  the  ray  optics  expression  eq.  (4.14)  virtually 
reproduces  the  curve  obtained  using  numerical  integration. 
If  one  uses  eq.  (4.18)  to  estimate  neff,  one  obtains  K values 

with  an  error  of  about  1%.  It  is  interesting  to  note  that 
the  linear  dependence  of  K on  ka  is  explicit  in  eq.  (4.14)  . 

In  the  standard  numerical  approach,  one  must  do  a number  of 
integrations  and  plot  the  data  properly  before  the  linear 
dependence  becomes  obvious . The  second  example  Figure  4-7  is 
a plot  of  K versus  the  thickness  of  the  active  layer  d for 

the  4-layer  structure  shown  by  Itaya  et  al.  [Ita84] . 
Excellent  agreement  between  the  two  approaches  is  observed 
over  the  whole  range  of  interest.  Using  eq.  (4.18)  to 
estimate  neff  for  eq.  (4.14),  one  gets  better  than  1%  accuracy 
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for  K when  t<0 . 2(im  and  d<0.3|lm.  For  t=0 . 3|lm,  the  error 
increases  to  about  5%  at  d=0.2|lm,  and  is  even  larger  for 
d>0.2p.m.  However,  the  region  where  d is  large  is  usually  of 
little  interest  to  DFB  laser  designers  since  threshold 
current  increases  as  d increases.  The  roll-over  behavior  for 
t— 0 . l|im  has  to  do  with  the  fact  that  B has  a maximum  at  some 


due  to  the  monotonic  decrease  in  0ef f • Since  K is  given  by  B 
, there  will  be  a peak  in  K near  this  d value.  Similar 

phenomena  also  occur  in  metalized  DFB  laser  configurations  in 
Figure  4-3.  In  the  next  example  Figure  4-8,  K is  plotted 

versus  the  number  of  quantum  wells  for  a 1.55jlm  MQW  DFB  laser 
as  shown  by  Makino  and  Adams  [Mak91] . The  ray  optics, 
coupled  mode  and  transfer  matrix  approaches  show  excellent 
agreement  for  such  a weak  perturbation.  In  this  MQW 
structure,  eq.  (4.18)  can  be  used  for  neff  in  the  ray  optics 
expression  without  introducing  any  noticeable  error  (much 
less  than  1%)  . It  is  interesting  to  note  that  there  is  a 
value  for  Nw  which  maximizes  K (for  the  present  case,  Nw=2) . 

This  is  also  similar  to  the  roll-over  behavior  in  Figure  4-7 . 
As  the  final  example,  Figure  4-9  is  a plot  for  K in  a DBR 

type  structure  with  a corrugated  semiconductor/air  interface 
where  the  perturbation  is  stronger.  The  ray  optics  formula 
agrees  with  the  transfer  matrix  method  very  good  in  this 
case.  Note  that  the  superlinear  behavior  is  because  the  t is 
assumed  to  be  reduced  when  the  grating  is  deeper.  The 


value  of 


monotonically  with  decreasing  d 


coupled  mode  results  are  inconsistently  underestimated  by 
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about  20-30%  at  only  ka=0.2,  depending  on  the  choice  of  the 
unperturbed  structure:  the  lower  dotted  curve  is  based  on  a 
4-layer  unperturbed  waveguide  where  the  nb/nc  interface  is  in 
the  center,  and  the  upper  dash-dot  curve  is  based  on  a 5- 
layer  unperturbed  waveguide  where  an  additional  "grating 
layer"  is  assumed  with  an  averaged  refractive  index. 

4 . 4 Mixed  Gain/Index  Coupled  Gratings 

So  far  the  coupling  coefficient  has  been  assumed  to  be 
real  and  only  the  magnitude  of  K has  been  calculated.  In  the 

most  general  case,  the  cladding  layer  may  be  lossy  due  to 
free  carrier  or  band-to-band  absorption  when  a metal  or  a 
small  bandgap  material  is  used.  In  such  cases  K becomes 
complex.  A complex  K can  in  principle  improve  the  single 
frequency  yield  [Kap82].  The  complex  coupling  coefficient 
should  be  written  as 

K = Ki  + iKg  = | K | (4.19) 
where  Ki  and  Kg  denote  the  index  and  gain  (or  loss)  coupling. 
For  practical  cases  with  loss  gratings,  Ki>>Kg  and  Ki= | K | 
since  one  does  not  want  too  much  loss  which  increases  the 
threshold  current.  Under  this  condition,  the  analysis 
presented  in  previous  sections  can  be  used  to  calculate  |K|, 
and  Kg  can  be  computed  in  two  ways. 

The  first  method  is  to  obtain  |K|  from  eq.  (4.12)  using 
proper  complex  refractive  indices.  Then  the  magnitude  of  Kg 
relative  to  Ki  is  determined  from  the  ratio  of  the  imaginary 
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O 

part  to  the  real  part  of  the  complex  index  perturbation  An 


for  the  grating  in  accordance  with  the  coupled  mode 
definition  eq.  (4.1). 

The  second  method  is  to  define  a complex  amplitude 
diffraction  efficiency  similar  to  eq.  (4.9)  but  without  the 
absolute  value.  Then  complex  K is  expressed  as 


Since  eq.  (4.20)  always  gives  a complex  number  even  for  a 
pure  index  coupling  structure,  one  has  to  find  the  reference 
phase  before  determining  the  actual  phase  <|)  of  the  complex 

coupling  coefficient.  Assuming  the  loss  is  small,  one  first 
turns  off  the  loss  term  by  setting  the  extinction  index  to 
zero  (for  dielectric)  or  infinity  (for  metal)  and  gets  Ko 
K0  = I Ko  I exp  < i<J)o ) 

Then  one  turns  on  the  loss  term  and  gets  Ki 
Ki  = | Ki  | exp  (i<(>i) 

Finally  the  complex  coupling  coefficient  is 

K = | Ki  | exp  ( i<J>i— i<>o ) (4.21) 

A properly  controlled  loss  coupling  DFB  mechanism  has 
been  shown  to  improve  the  single  frequency  yield  [Tsa92] . In 
the  shiny  contact  DFB  laser,  the  loss  coupling  can  be 
controlled  by  the  metalization  process  since  the  shiny  metal 
(gold)  is  almost  lossless.  For  example,  thin  layers  of 


K 


] 


(4.20) 


2 neff  weff 


titanium  and  platinum  can  be  deposited  first  on  the 
corrugated  semiconductor  for  better  p-contact  reliability. 
Since  titanium  is  lossy  and  not  a shiny  metal,  this  layer 
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gives  rise  to  periodic  loss  modulation.  Alternatively,  a 
thin  contact  layer  of  heavy  doping  may  be  grown  on  top  of  the 
cladding  layer  for  better  electrical  resistance.  The  free 
carrier  or  band-to-band  absorption  in  this  contact  layer 
introduces  a loss  coupling  DFB  mechanism.  Optimization  of 
these  layer  thicknesses,  which  can  be  done  using  the  extended 
ray  optics  approach,  should  lead  to  DFB  lasers  of  high 
reliability  and  high  single  frequency  yield. 
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Figure  4-1:  Diagram  for  the  ray  optics  calculations  of  K 
for  a generic  DFB  laser  structure 


coupling  coefficient  k (cm-1) 
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normalized  groove  depth  ka 


Figure  4-2:  The  coupling  coefficient  for  the  DFB  laser 
in  Figure  3-1  with  the  following  parameter  values: 

A.=  0 . 85|lm,  d=0  . l|lm,  t=0 . 3|lm,  na=3 . 6,  nb=3 . 4 , nd=3 . 4 


coupling  coefficient  k (cm-1) 
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buffer  layer  thickness  t (pm) 


Figure  4-3:  The  coupling  coefficient  for  the  DFB  laser 
in  Figure  3-1  with  the  following  parameter  values: 
^=0 . 85|lm,  d=0  . l|lm,  ka=0 . 2,  na=3 . 6,  nb=3 . 4,  nd=3 . 4 


coupling  coefficient  k (cm-1) 


68 


active  layer  thickness  d (|im) 


Figure  4-4:  The  coupling  coefficient  for  the  DFB  laser 
in  Figure  3-1  with  the  following  parameter  values: 
X=0 . 85^lm,  ka=0 . 2,  na=3 . 6,  nb=3 . 4 , nd=3 . 4 
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Figure  4-5:  Generic  DFB  structure  with  dielectric  gratings 

(a)  SCH  DFB  laser  configuration 

(b)  MQW  DFB  laser  configuration 
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normalized  groove  depth  ka 


Figure  4-6:  The  coupling  coefficient  for  the  SCH  DFB  laser 
in  Figure  4-5  (a)  for  sinusoidal  gratings 
with  the  following  parameter  values: 
na=3 . 54 , nb=nb '=3.40, nc=nc '=3.17 
A,=  l . 55|lm,  d=0  . l|lm,  t=t ' =0  . l|im 


coupling  coefficient  k (cnr1) 
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active  layer  thickness  d (jam) 


Figure  4-7:  The  coupling  coefficient  for  the  SCH  DFB  laser 
in  Figure  4-5  (a)  for  sinusoidal  gratings 
with  the  following  parameter  values: 
na=3 . 54 , nb=3 .39, nc=nc '=3.17 
A,=  l . 55(lm,  t ' = 0, 2a=500A 


coupling  coefficient  k (cm'1) 
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number  of  quantum  wells 


Figure  4-8:  The  coupling  coefficient  for  the  MQW  DFB  laser 
in  Figure  4-5 (b)  for  square  gratings 
with  the  following  parameter  values: 
nqW— 3 . 57 , nbr=3 . 39, nb=nb '=3 . 39, nc=nc '=3.17 
dqw=70A,  dbr=100A,  t = 0 . 185|lm,  t ' =0  . Ol^lm 
A,=  l . 55|lm,  2a=300A 


coupling  coefficient  k (cm-1) 
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normalized  groove  depth  ka 


Figure  4-9:  The  coupling  coefficient  for  a DBR  structure 
with  square  gratings  and  the  following  parameter  values: 
X=1 . 8 Jim,  d=  1 Jim , t+a=l|lm,  t ' =0 , na=5 . 0 , nb=4 . 7 , nc=l . 0,  nc  ' =4 . 7 


CHAPTER  V 

DEEP  GROOVE  GRATINGS 


In  previous  chapters,  only  shallow  groove  gratings 
(ka«l)  have  been  considered  in  the  analysis.  Such  gratings 
are  commonly  utilized  in  DFB  lasers  in  fiber  communication 
applications.  In  this  chapter,  we  discuss  the  theoretical 
aspects  of  deep  groove  gratings  and  propose  a grating  surface 
emitter  (GSE)  structure  with  a deep  groove  metalized  grating. 
This  shiny  contact  GSE  structure  can  be  potentially 
fabricated  into  two  dimensional  arrays  of  low  threshold 
current,  single  frequency  devices  competitive  with  vertical 
cavity  surface  emitting  lasers  (VCSEL)  [Jew91] . In  Section 
5.1,  the  formalism  for  the  diffraction  efficiency  for  deep 
groove  gratings  is  outlined  using  the  Rayleigh  approach. 
Then  numerical  results  for  first  and  second  order  gratings 
are  given  in  Sections  5.2  and  5.3,  and  the  design  of  such 
potential  GSE  devices  discussed. 

5 . 1 General  Formalism 

In  this  section,  the  shiny  contact  DFB  laser  as  shown  in 
Figure  5-1  is  discussed.  This  DFB  structure  is  similar  to 
the  one  drawn  in  Figure  3-1  except  for  arbitrary  groove 
depth.  Assume  the  cladding  layer  is  made  of  sinusoidally 
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corrugated  perfect  metal  and  we  consider  TE  modes,  i.e.,  the 


definition  for  the  coupling  coefficient  K in  Chapter  IV, 


Now  we  show  how  to  compute  the  diffraction  efficiency  T)  for  a 
grating  with  depth  comparable  to  or  even  larger  than  the 
wavelength  (ka~l) . As  mentioned  in  Chapter  III,  even  though 
the  Rayleigh  assumption  is  questionable  for  the  general  case, 
the  approach  itself  is  straightforward  and  can  provide  fast 
and  accurate  results  for  sinusoidal  gratings  compared  with 
other  techniques  [Pet80]  . Therefore  we  use  the  Rayleigh 
expansion  in  this  chapter. 

Assume  the  incident  wave  has  unit  amplitude  and  x-  and 
z-propagation  constant  (Jo  and  po  respectively  (Po=neffk)  : 

exp  (iOox)  exp  (~i(30z) 

Then  the  total  electric  field  in  layers  na  and  n^  can  be 
expressed  as,  according  to  Rayleigh  assumption. 


electric  field  is  y-oriented.  Recall  the  ray  optics 


(5.1) 


OO 


Ea(x,z)  = exp  (iOox)  exp  (-i|30z)  + ^ Amexp  (~iamx)  exp  (-ifimz) 


m=-oo 


OO 


Ec(x,z)  = 0 in  the  perfect  metal 


(5.2) 


Pm  = Po  + mK  = p0  + m(2jt/A)  A:  grating  periodicity 


where 
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Pm  = VPm2-nb2k2 

k = 27C/A,  X:  wavelength  in  vacuum 

Care  should  be  taken  for  the  sign  of  am . The  former 
expression  is  used  for  na>neff  and  the  latter  for  na<neff. 
Following  a boundary  matching  procedure  similar  to  that 
presented  in  Section  3.1,  the  infinite  inhomogeneous  linear 
system  is  so  obtained, 

Aoexp  (-iaot ) - Boexp(-pot)  - Coexp(pot)  = -exp(iaot) 
idoAoexp (-iGot)  - poBoexp (-pot)  + poCoexp(pot)  = iOoexp(iOot) 
Amexp (-iomt)  - Bmexp(-pmt)  - Cmexp(pmt)  = 0 Vm*0 

iCTjnAmexp  ( — iCJmt)  - Pi^Bj^exp  ( — Pmt ) + PmCmOxp  (Pmt ) = 0 Vm^O 

oo 

[ (~1)  BmIm_n  (pma ) + CmIm-n(Pma)  ] = 0 Vn 

m=-oo 

(5.3) 

where  In  is  the  modified  (or  hyperbolic)  Bessel  function  of 
the  first  kind  (of  order  n) . This  system  can  be  solved  using 
standard  numerical  techniques,  and  the  results  for  two 
special  cases  are  given  below. 

5 . 2 First  Order  Deep  Groove  Gratings 

For  a first  order  grating,  K=2(3o  or  A=A./ (2neff)  , there 
are  two  traveling  terms  that  carry  electromagnetic  power:  the 
fundamental  forward  wave  (m=0,  amplitude  Ao , propagation 
constant  Po)  and  the  backward  diffracted  wave  (m=-l, 
amplitude  A_i,  propagation  constant  -(3q)  . The  diffraction 


77 


efficiency  is  then  T|=|A_il2,  and  numerical  data  for  the 
coupling  coefficient  K can  be  obtained  from  eq.  (5.1)  . 

In  order  to  see  how  many  terms  are  required  for  the 
Rayleigh  expansion  eq.  (5.2)  to  converge,  K is  plotted  in 

Figure  5-2  versus  the  number  of  terms  in  the  expansion  for 
ka=0.5  using  the  parameter  set  as  in  Figure  4-2.  From  this 
figure  it  is  obvious  that  about  10  terms  are  needed  for  a 
stable  and  convergent  output,  and  this  corresponds  to  a 30x30 
matrix.  Then  we  plot  the  results  for  K versus  ka  in  Figure 

5-3  for  the  same  structure,  where  the  solid  curve  is  obtained 
using  the  numerical  data  for  T]  from  an  expansion  of  at  least 

10  terms,  dashed  line  from  the  closed-form  expression  eq. 
(4.13)  based  on  a two  wave  approximation,  and  the  dotted  line 
from  the  transfer  matrix  method.  It  is  observed  that  for 
small  values  of  ka,  there  is  no  significant  difference  among 
these  three  curves.  For  larger  values  of  ka,  however,  the 
accurate  results  should  be  sublinear  and  then  rolling  over. 
The  two  wave  approximated  eq.  (4.13)  produces  linear  results 
while  the  transfer  matrix  goes  super  linear.  This  indicates 
that  the  ray  optics  formalism  is  superior:  for  best  accuracy, 
one  can  use  numerical  data  for  T|,  and  for  simplicity,  one  can 

use  the  closed-form  expression  which  is  still  better  than  the 
other  techniques. 

The  results  in  this  section  show  that  for  a first  order 
grating,  K does  not  always  go  up  when  the  groove  depth 
increases.  Therefore  one  has  to  be  careful  in  maximizing  K 

using  a deep  groove  device.  There  is  an  optimal  value  of  the 
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groove  depth  for  maximum  diffraction  efficiency.  Similar 
behavior  for  the  diffraction  efficiency  of  a corrugated 
air/metal  grating  has  been  reported  in  the  literature  using 
other  techniques  such  as  [Moh86] . 

5 . 3 Second  Order  Deep  Groove  Gratings 

For  a second  order  grating,  K=Po  or  A=A./neff,  the  basic 

three  terms  are  m=0,  m=-l  and  m=-2  whose  z-propagat  ion 
constants  are  Po,  0 and  -po  respectively.  The  first  order 

diffracted  wave  A_i  corresponds  to  the  surface  emitting  beam 
and  the  second  order  diffracted  wave  A_2  provides  the  backward 
coupling . In  a low  threshold  GSE  design  when  low  surface 
emitting  power  is  needed,  one  wants  to  maximize  A_2  but 
minimize  A_i  since  A_2  provides  the  feedback  but  A_i  is  the 
radiation  loss  to  the  laser  cavity.  Such  a design  can  not  be 
realized  using  simplified  theories  valid  only  for  shallow 
groove  gratings.  For  example,  when  higher  order  waves  are 
neglected,  A_2  is  always  much  less  than  A_i  [Rig76], 
consequently  conventional  GSE  lasers  using  shallow  groove 
gratings  show  high  threshold  current  while  produce  high 
surface  emitting  power,  and  this  is  why  they  have  not  been 
utilized  in  low  threshold,  low  power  applications  such  as 
chip-to-chip  communications. 

In  this  section  we  investigate  deep  groove  effects  in 
second  order  metalized  gratings  by  including  those  higher 
order  evanescent  waves  in  the  ray  optics  formalism.  The 
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results  presented  in  this  section  will  shed  light  on 
potential  low  threshold  GSE  lasers  which  could  be  competitive 
with  VCSELs  [Tak92 ] . 

Because  the  grating  wavevector  is  smaller  in  the  second 
order  grating,  higher  order  waves  do  not  decay  as  fast  as  in 
the  first  order  case.  We  find  that  about  15  terms  are 
required  for  convergent  results,  and  this  corresponds  to  a 
45x45  matrix.  Both  |A_il  and  |A_2l  are  plotted  in  Figure  5-4 

for  the  same  configuration  as  in  Figure  5-3  except  for  the 
second  order  grating.  In  the  region  ka<0.5,  | A_i | is  much 
larger  than  | A_2  I as  expected,  and  a large  portion  of  the 
optical  power  is  emitted  as  surface  radiation  loss.  For  high 
power  GSE  lasers  where  the  threshold  is  less  concerned,  a 
large  |A_il  is  desired  for  large  surface  emitting  power.  For 
low  threshold  GSE  or  edge  emitting  DFB  lasers,  however,  a 
large  I A_  2 I for  the  feedback  and  a small  |A_il  for  low 
radiation  loss  are  preferred.  This  cannot  be  achieved  using 
shallow  groove  gratings,  so  most  people  do  not  use  second 
order  gratings  for  edge  emitting  DFB  lasers  even  though  the 
single  frequency  yield  is  better  [Car92].  When  the  groove 
depth  becomes  bigger,  |A_il  starts  to  roll  over  while  | A_2 1 
keeps  rising.  At  around  ka=l,  |A_il  is  much  less  than  | A-2 1 , 
and  this  is  the  region  of  interest  for  such  a novel  low 
threshold  GSE  or  DFB  laser.  For  a low  power  GSE  laser,  one 
can  in  principle  achieve  threshold  currents  as  low  as  for 
conventional  DFB  lasers.  For  a DFB  laser  using  second  order 
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deep  groove  gratings,  one  should  be  able  to  demonstrate  high 
single  frequency  yield  and  easy  fabrication. 

The  oscillatory  behavior  in  Figure  5-4  using  the 
Rayleigh  expansion  shows  agreement  with  the  diffraction  of  a 
corrugated  air/metal  interface  using  other  techniques 
[Moh86]  . We  have  also  observed  in  our  laboratory  a similar 
phenomenon  that  high  order  waves  are  brighter  than  the  low 
order  ones  when  the  grating  depth  is  comparable  to  or  larger 
than  the  wavelength.  Therefore,  it  is  believed  that  such  GSE 
lasers  using  deep  groove  shiny  contact  gratings  can  be 
realized  with  a proper  design.  Due  to  the  superior  thermal 
properties  inherent  in  such  deep  groove  shiny  contact  GSE 
laser  structures,  they  could  be  competitive  with  vertical 
cavity  surface  emitting  lasers  in  future  chip-to-chip 
communication  technologies. 
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Figure  5-1:  Diagram  for  the  ray  optics  calculations  of  K 
for  a generic  DFB  laser  structure 


coupling  coefficient  k (cm'1) 
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number  of  terms  in  Rayleigh  expansion 


Figure  5-2:  The  coupling  coefficient  for  the  DFB  structure 
in  Figure  5-1  and  the  following  parameter  values: 

A.=  0.85|lm,  d=0.1(im,  t=0.3flm,  ka=0.5 
na=3.6,  nb=3.4,  nc=-i<*>,  nd=3 . 4 


coupling  coefficient  k (cm-1) 
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normalized  groove  depth  ka 


Figure  5-3:  The  coupling  coefficient  for  the  DFB  structure 
in  Figure  5-1  and  the  following  parameter  values: 

A.=  0.85(lm,  d=0.1flm,  t = 0.3|lm, 
na=3.6,  nb=3.4,  nc=-i°°,  nd=3 . 4 


relative  amplitude 
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normalized  groove  depth  ka 


Figure  5-4:  The  relative  amplitude  for  two  diffraction  orders 
of  a second  order  DFB  structure  with  these  parameter  values: 
X=0 . 85|lm,  d=0  . l|lm,  t=0.3|im,  na=3.6,  nb-3.4,  nd=3 . 4 


CHAPTER  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


In  conclusion,  shiny  contact  diode  lasers  have  been 
demonstrated  for  the  first  time  in  the  1 . 3(lm  InGaAsP/InP 
material  system.  Power  levels  greater  than  one  watt  pulsed 
have  been  achieved  at  room  temperature,  among  the  highest  for 
DH  InGaAsP/InP  lasers  and  comparable  to  those  obtained  from 
state-of-the-art  separate  confinement  heterostructure  single 
quantum  well  lasers.  The  experimental  data  have  been 
successfully  explained  by  the  theory  and  key  parameters 
determined  from  the  measurements.  Using  measured  parameters, 
improvements  in  continuous  wave  performance  levels  relative 
to  conventional  devices  are  predicted.  We  believe  that  these 
improvements  can  be  realized  according  to  our  optimal  design, 
and  it  is  recommended  for  future  research. 

The  advantages  of  combining  the  shiny  contact  with  the 
DFB  laser  structure  have  been  made  clear.  Ray  optics  and 
Floquet-Bloch  formalisms  have  been  developed  to  model  the 
lasing  properties  of  these  shiny  contact  DFB  lasers  as  well 
as  conventional  DFB  lasers  . It  has  been  shown  that  in 
standard  type  DFB  lasers  with  dielectric  gratings  where  the 
index  difference  is  small,  the  K values  from  the  coupled  mode 

overlap  integral,  the  transfer  matrix  method  and  the  ray 
optics  expression  agree  very  well.  In  shiny  contact  DFB 
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lasers  with  metalized  gratings  where  the  index  difference  is 
very  large,  the  standard  techniques  do  not  work  well  and  the 
ray  optics  expression  is  consistent  with  the  Floquet-Bloch 
approach.  The  ray  optics  formalism  has  been  extended  to 
investigate  deep  groove  effects  of  second  order  metalized 
gratings.  Our  calculations  indicate  that  low  threshold 
grating  surface  emitting  laser  arrays  are  possible  by  the  use 
of  deep  groove  shiny  contact  gratings.  GSE  arrays  of  this 
type  could  be  competitive  with  vertical  cavity  surface 
emitting  laser  arrays  due  to  their  superior  thermal 
properties . 

For  future  research,  the  shiny  contact  DFB  laser 
structure  with  the  loss  coupling  mechanism  mentioned  in 
Chapter  IV  offers  several  advantages  absent  in  conventional 
DFB  lasers.  We  believe  that  it  is  worthwhile  to  make  such 
devices  in  the  InGaAsP/InP  system  since  the  key  technology, 
fabrication  of  the  shiny  contact,  has  been  realized  in  this 
dissertation.  The  low  threshold  GSE  laser  array  using  deep 
groove  shiny  contact  gratings,  first  proposed  in  this  study, 
is  also  very  interesting.  Experimental  verification  of  such 
devices  could  have  a potential  impact  on  future  chip-to-chip 
communication  technologies. 
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